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Radioactive Decay 


DECUS Program Library Write-up 


DECUS NO. 8-542 


1) Program 

This program solves the radioactive decay equation; 

A=Ao , A=(A0)*EXP(-DT) 

for any one of the four variables: 

A=activity after decay 

A0=activity before decay=original activity 
TH=half-life of the radioisotope 
T=time of decay,where: 

D=the decay constant(A)= .69515/TH. 


2) Operation 

The program is written in FORTRAN-D and is loaded by 
FOSL.See DEC'S Introduction to Programming(1970) p.8--40. 

When started the program requests the code of the 
variable to be calculated.The code is as follows; 

A=1 A0=2 T=3 TH=4 

The program then requests the number of calculations 
to be made,A maximum of 30 is permitted. 

Depending on the variable to be calculated the program 
then requests the values of the other three variables to be 
inputted. 

When the input data is complete the program outputs, 
in tabulated form,the input data together with the calculated 
variable. 

The program is endless i.e. after completing a set of 
computations it re-starts itself. 











3) Sample Run 


.FOSL 

♦IM-R: 

* 

+OPT- 

t 

♦READY 

t 


TYPE CODE OF VARIABLE TO BE CALCULATED 
A=1 A0=2 T=3 TH=A : 1,,__— 


typed 1 followed by a space 


HOW MAHY CALCULATIOMS 12 . ___ tvned 2 followed by a space 

TYPE A0>T>TH ^ 

10000. 8.02 8.02 ) typed a value for A0, space,value for T,space, 

100. 8.02 8.02 ) value for TH and carriage return 


TH 

0.802000E+1 

0.802000E+1 


T 

0.802000E+1 
0.802000E+1 


A0 

0.100000E+5 
0. 100000E+3 


A 

0.A99999E+A 
0.A99999E+2 


TYPE CODE OF VARIABLE TO BE CALCULATED 
A=1 A0=2 T=3 TH=4 :2 

HOW MANY CALCULATIONS ?2 
TYPE A>T#TH 
5000 8.02 8.02 

50 8.02 8.02 


TH 

0.802000E+1 

0.802000E+1 


T 

0.802000E+1 

0.802000E+1 


A0 

0.100000E+5 
0.100000E+3 


A 

0.500000E+4 
0.500000E+2 


TYPE CODE OF VARIABLE TO BE CALCULATED 
/V=l A0=2 T=3 TH=4 :3 

HOW MANY CALCULATIONS ?2» 

TYPE A>A0#TH 
100 200 8.02 

5000 10000 8.02 


TH 

0.802000E+1 

0.802000E+1 


T 

0.801997E+1 

0.801995E+1 


A0 

0.a00000E+3 
0. 100000E+5 


A 

0.100000E+3 
0.500000E+4 


TYPE CODE OF VARIABLE TO BE CALCULATED 
A=1 A0=2 T=3 TH=4 : 












4) Characteristics 


a. All input and output is via the teletype. 

b. The units of time T and TH must be identical i.e. seconds, 
minutes,days etc... 

c. The units of activity in A and A0 must be identical i.e. 
Ci(Curie),uCi,counts/sec. etc... 

d. The variables code (lto4) and the number of calculations 
(up to 30) are integers.However,they can be terminated 
by a space,comma,period,carriage return etc... 

e. For clarity it is best to terminate each group of three 
input data by a carriage return and terminate an individual 
number by one or more spaces. 

f. Core storage for the program and the FOETEAN-D operating 
system is 4K. 









5) Listing 


C 

C 

c 

c 

c 

c 

c 

c 


1 


30 

32 


35 


38 


RADIOACTI 
TO SOLVE 
FOR A>A0> 
A =ACTIVI 
A0=ACTIVI 


VE DECAY BY A* L•AL-MOAIMI»APRIL 1972* 
THE DECAY EQUATION A=(A0)*E<P(-DT) 

T OR TH WHERE* 

TY AFTER DECAY. 

TY BEFORE DECAY=ORI 61MAL ACTIVITY. 


D “DECAY COMSTAMT=.69315/TH. 
TH=HALF LIFE OF THE ISOTOPE. 
T “TIME OF DECAY. 


DIMEMSIOM TH< 30)> T( 30)# A0(30)> A(30)> G(30)> E(30) 

TYPE 6 
TYPE 7 
ACCEPT 8>L 
TYPE 20 
ACCEPT 


GO TO <30»A0#50#60)>L 
TYPE 9 
DO 32 I“1#M 

ACCEPT 1 5j A0(I)>T<I)>THII) 

COMTIMUE 

DO 35 I“1>M 

G( I ) = . 6931 5*T( I )/TH< I ) 

EC I )=EXPF( -G( I ) ) 
A(I)“A0<I)*ECI) 

CONTIMUE 
TYPE 16 
DO 38 I“1*N 

TYPE 17>TH(I)>TCI)>A0CI)»ACI) 

CONTINUE 

GO TO 1 


40 TYPE 10 

DO 42 I“1>N 

ACCEPT 15jAC I)>T(I)*THCI) 

42 CONTINUE 

DO 45 I“UN 

GC I ) “ . 6931 5i«TC I ) /THC I ) 

EC I ) “EXPFC -GC I ) ) 

A0CI)=ACI)/EC I) 

45 CONTINUE 

TYPE 16 
DO 48 I“1>N 

TYPE 17»THCI)#TCI)>A0CI).AC I) 
48 CONTI N.UE 

GO TO 1 


4 









5) Listing (continued) 


50 TVPE 11 

DO 52 I=1>M 

ACCEPT 15*AC IA0<I)>THCI) 

52 CO!>JTI'MUE 

DO 55 I = 1>^J 

G(I)=.69315/TH<I) 

EC I ) =LOGFC A0C I) ) -LOGFC AC I ) ) 

TCI)=ECI)/GCI) 

55 COMTIMUE 

TTPE 16 
DO 58 I = 1#M 

TYPE 17»THCI)>TCI)»A0CI>>AC I) 

58 COMTIMUE 

GO TO 1 

60 TYPE 12 

DO 62 1=1#M 

ACCEPT 15 #AC I)»A0C I)#TCn 
62 COMTIMUE 

DO 65 I=1#M 
GC I) = .69315*TC I) 

EC I ) =LOGFC A0C I ) )-LOGFC AC I ) ) 

THCI)=GCI)/ECI) 

65 CONTINUE 

TYPE 16 
DO 68 I = 1#M 

TYPE 17# THCI)# TCI)#A0CI)#AC I) 

68 CONTINUE 

GO TO 1 

6 FORMATC/////"TYPE CODE OF VARIABLE TO BE CALCULATED") 

7 FORMATC/"A=l A0=2 T=3 TH=A :”) 

8 FORMATC I ) 


9 

FORMATC/"TYPE 

A0#T#TH"/) 

10 

FORMATC/"TYPE 

A# T# TH"/) 

11 

FORMATC/"TYPE 

A# A0#TH"/) 

12 

FORMATC/"TYPE 

A#A0#T"/) 

15 

FORMATC E#E#E) 


16 

FORMATC/" 

TH T 

17 

FORMATC/E# E# E# 

E) 

20 

FORMATC/"HOW MA^MY CALCULATIONS ?") 


END 



♦ 
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